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Location, correlation, radiation: 
where is the cr, what is its structure and what is its coupling to photons? 
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Scalar mesons are a key expression of the infrared regime of QCD. The lightest of these is the cr. Now 
that its pole in the complex energy plane has been precisely located, we can ask whether this state is 
transiently qq or qqqq or a multi-meson molecule or largely glue? The two photon decay of the cr can, 
in principle, discriminate between these possibilities. We review here how the yy — > n*7t" , n°n° cross- 
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sections can be accurately computed. The result not only agrees with experiment, but definitively fixes 
the radiative coupling of the cr. This equates to a two photon width of (4.1 ± 0.3) keV, which accords 
with the simple non-relativistic quark model expectation for a uu , dd scalar. Nevertheless, robust 
predictions from relativistic strong coupling QCD are required for each of the possible compositions 
before we can be sure which one really delivers the determined yy coupling. 

Keywords: Scalars; photon couplings; QCD 

1. Location: the vacuum of QCD and the spectrum of scalars 

tur 
1 



The nature and role of scalar mesons is one of the most intriguing enigmas of strong dy- 
namics - . The scalars are special: they constitute the Higgs sector of the strong interaction, 
giving masses to all of the light hadrons. They are in turn intimately tied to the structure of 
5_i ■ the QCD vacuum. To understand this we have to go back 60 years. We know that just count- 

ing quarks is sufficient to explain why a typical meson, like the p, and the lightest baryon, 
the nucleon, have masses in the ratio of 2 : 3. Then the question is why are pions so very 
light? Indeed, since the current masses of the up and down quarks are just a few MeV, one 
could equally ask why are the masses of typical mesons and baryons so heavy! The answer 
is in the structure of the QCD vacuum. The vacuum of QED is essentially empty with just 
perturbatively calculable corrections from virtual photon emission and absorption, includ- 
ing e + e~ pairs. The vacuum of QCD is quite different. While asymptotic freedom ensures 
that quarks propagate freely over very short distances deep inside a hadron, little perturbed 
by the sea of quark-antiquark pairs and the clouds of gluons, over larger distances (distances 
of the order of a fermi), the interactions between quarks and gluons become so strong that 
they produce correlations in the vacuum. Correlations that we know as condensates. It is 
by moving through this complex vacuum that the light quarks gain their mass. 

We are, of course, familiar with the idea that the world of hadron physics reflects the 
symmetry properties of the underlying world of quarks and gluons. The hadron world has 
isospin symmetry because of the near equal mass of the up and down quarks and their 
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flavour blind interactions. The fact that the current masses of these quarks are so much 
lighter than the natural scale of QCD, viz. Aqcd, means that we can regard them as mass- 
less. Then their helicity becomes a well-defined quantum number and the left and right- 
handed worlds decouple. QCD has an almost exact S U(2) x S U(2) chiral symmetry. In 
the world of hadrons this symmetry is not apparent. Scalars and pseudoscalars, vectors and 
axial-vectors are not degenerate in mass with closely related interactions. Consequently this 
symmetry must be broken. Long long ago before QCD was uncovered, Nambu proposed a 
mechanism for this spontaneous breakdown of chiral symmetry I ^PH I jf me potential gen- 
erated by the interactions of a scalar (<x) field and the pseudoscalar (n) has a Mexican hat 
shape, then nature chooses a ground state in which pions correspond to the quantum fluc- 
tuations round the bottom of the hat and are massless. In contrast, the physical scalar field, 
produced by fluctuations about its non-zero vacuum expectation value are up and down the 
sides of the hat, implying it is massive. What is happening in the underlying world of QCD 
is that the operators qq, qGq, etc., gain non-zero vacuum expectation values. To realise 
Nambu 's picture, this dynamical breaking must be dominated by the qq condensate. 

This is all very well-known. What is relatively new, is that these ideas have been tested 
in experiment and found to agree exactly with this picture. For twenty years QCD sum- 
rules have indicated that the value of (qq) -(240 MeV) 3 and it is this 240 MeV that 
sets the scale for the constituent mass of the up and down quarks. Similarly, Schwinger- 
Dyson calculations of the Euclidean momentum dependence of the quark mass function 
show that such a value for the chiral condensate follows just from Aqcd- Because pions 
are the Goldstone bosons of chiral symmtery breaking, their interactions at low energy 
directly reflect the value of this condensate. Such low energy interactions can be studied 
experimentally in K e \ decay for instance. There a K decays into an electron and its neutrino, 
plus two pions, which once formed interact strongly. Because pions are the lightest of 
all hadrons, their interactions are universal. Thus from the K e \ decay distribution we can 
learn about the phase of the relevant nn interaction. The decay depends on five kinematic 
variables El the di-lepton mass, the di-pion mass, the three angles, Of, 9 n and <p. In the 
kaon rest-frame the direction of the back-to-back di-lepton and di-pion systems defines 
a direction. Then in the di-lepton rest frame, each lepton goes off at an angle, 9(, to this 
direction, and the pions in the di-pion rest frame at an angle 9 n . The plane of the two 
leptons and the plane of the two pions are at angle (f> to each other. The 5-dimensional 
decay distribution then determines the relative phase of the nn S and P-wave interactions 
as a function of di-pion mass. The first high statistics measurement came 6 or 7 years ago 

from E865 at BNL si These showed that the low energy behaviour of this phase difference 

— l7l — 

corresponded to a qq condensate of exactly the size expected 1 -^. This dominance of the qq 

condensate has been further checked by the preliminary results from NA48^1 

Thus we see that Nambu's picture of the breaking of chiral symmetry is indeed realised 

in nature. So what is the scalar field, the non-zero vacuum expectation value of which 

generates masses for all light hadrons? Do we identify Nambu's <x with the lightest / = 

7 = state, the /o(600) ? That such a meson might exist was proposed more than 40 

years ago. However, Breit-Wigner fits to nn mass distributions have found a wide range of 

masses and widths over the years, making it far from certain that such a state really existed. 
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Fig. 1. The complex energy squared (s) plane is displayed. The position of the T-matrix pole for the <x deter- 
mined in each of the analyses listed in the P PG table s is show n. The position of the p-pole is marked for 
orientation. Singled out are those labelled ProtDS BFpGS EstGSand AMpHlfrom 30 year s ago. CCL from the 
more recent Roy equation analysis and BES from a fit to the BES J/i// — » a>{nit) data — are highlighted. The 
PDG listings IS1 should be seen for the complete list of other references. 



For instance, the AMP analysis , the most complete in the mid '80's, found a pole at 
E = 0.87 - /0.37 GeV. A later treatment by Anisovich and SarantsevES which included fits 
to the production of two pion subsystems in ~pp annihilation from Crystal Barrel, found no 
pole whatsover. In contrast, more recently simple Breit-Wigner fits to n + n~ final states in 
D + — > n + n~n + data from E791 ^ and in J /iff — > u>n + n~ from BESII 021 have found poles 
at E = (478 ± 29) - /(162 + 1 1) MeV and E = (541 + 29) - /(252 + 42) MeV, respectively. 
However, these fits for the most part ignore the related phase variation determined for 
nn — > nn in the elasti c region. The scatter of pole positions from all the analyses quoted in 
the latest PDG tables EEHMMn lEl 

are shown in the complex s — E 2 plane in Fig. 1 . 
Now we know, thanks to a rather precise analytic continuation into the complex energy 

rm 

plane by Caprini, Colangelo and Leutwyler LL -4 that a pole does exist, signalling a state in 

the spectrum of hadrons. The analytic continuation starts with fixed-f dispersion relations. 

As always the left hand cut requires knowledge of crossed-channel exchanges. The unique 

feature of nn scattering is its 3-channel crossing symmetry. This means both the right and 

left hand cuts involve integrals over nn — > nn scattering. Making use of this crossing 

1 R 

property, and then partial wave projecting, leads to the Roy equations °. Solving this set 
of rigorous equations inputting experimental information above 800 MeV, together with 



I 



I 



1,2008 8:4 WSPC/1NSTRUCTION FILE review-mpla-rev 



4 M.R. Pennington 



ImE 




Re E 



Fig. 2. An illustration of the sheets and cut structure of the complex energy plane in a world with just one 
threshold and how these are connected. This represents the structure relevant to nil scattering near its threshold. 
Experiment is performed on the top shaded sheet, just above the cut along the real energy axis. The cross on the 
lower, or second, sheet indicates where the cr-pole resides at E 2 = s = sr. The ellipse above this on the top, or 
first, sheet indicates where the 5 -matrix is zero. 

chiral constraints, allows the / = J — nn partial wave to be determined everywhere on 
the first sheet of the energy plane Fig. 2. As shown by Caprini et al. 021 this fixes a 
zero of the S -matrix (symbolically depicted by the solid ellipse in Fig. 2) at £ = 441 - 
i 227 Me V, which reflects a pole (denoted by the cross) on the second sheet at the same 
position. This not only confirms the cr as a state in the spectrum of hadrons but locates 
the position of its pole very precisely with errors of only tens of MeV. This is within the 
region found by Zhou et al , who also took into account crossing and the left-hand cut. 
While the chiral expansion of amplitudes can have no poles at any finite order, particular 

im 

summations may. The inverse amplitude method 1 ^ 1 is one such procedure. Application of 
this by Pelaez et al.'^'did indeed find a pole in the same domain as Caprini et al. but many 
years earlier. However, without a proof that the Inverse Amplitude Method, rather than 
any other, provided precision unitarisation of the low order chiral expansion, the present 
author rather believed the analysis of a wide range of data of Ref. 10 that indicated no 
pole (or perhaps a very distant one). Now we know differently. There is a pole with a 
we ll-defined loc ation. This is far from the position proposed by the treatment of Ishida et 
al. | 23 | 24 | 25 | 26 | 27 |^ me deficiencies of which were explained long ago in Ref. 21. (BuggES 
has added to these arguments in response to the discussion on the position of the k by the 
same groupESl.) 

With a narrow resonance, there would naturally be a close correlation between the phase 
variation of the underlying amplitude on the real axis and in the complex plane, as one 
passes the pole. However, for the very short-lived cr, "deep" in the complex energy plane, 
this simple connection is lost. In Fig. 3 we show the phase of the / = J — nn — > nn 
scattering amplitude along two lines in the complex energy plane. The phase along the 
real axis is compared with relevant data from scattering 122-1 and K e 4 decays^ in Fig. 3. 
One sees how different this phase is compared with that on the lower sheet of Fig. 2 at 
ImE = -0.25 GeV. That deep in the complex plane shows the 180° phase change expected 
of a resonance. It is the dramatic variation in the amplitude as one moves away from the 
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real axis that has made the cr's existence and location so uncertain for so long. With the 
right tools for accurate analytic continuation, this is no longer an issue. 

But is this cr the chiral partner of the n&b i s it the Higgs of the strong interaction? 
There exist many other isoscalar scalars / (980),/o(1370),/o(1510),/ (1720). What role 
do these play? What is more, these have isodoublet and isotriplet partners, like the ,7^(800) 
(or k), A"p(1430), «o(980) and a ( )(1430). How are these all related? The quark model would 
lead us to expect a 3 Po nonet, but there are 19 scalars just listed: enough for two nonets 
with an isosinglet left over! Of course, there is a danger of subjectivity in what is included 
amongst this list. That the nine lightest do not fit a simple qq pattern anticipated from vector 
and tensor mesons is readily seen by noting that both the /o(980) and ao(980) sit very close 
to the KK threshold and couple very strongly to these channels. In a simple quark multiplet 
only the ss state does that. 

The possibility of discovery of pentaquark states has revived interest in c olour 3 scalar 
diquarks, formed from quarks of different flavour, like [ud], [us] and [ds] I34|35 [ Such 
diquarks can bind with the corresponding anti-diquarks to form tetraquark systems. The 
lightest multiplet is a nonet and lighterl^than the corresponding 3 Pq qq nonet, as in Fig. 4. 
Moreover, the heaviest tetraquarks are an /o and a^: [ns] [ns], with n — u,d, suggesting these 
might be good assignments for the 980 MeV states. However, Weinstein and i S g Ur l26| 37 | 38 | 
found within their potential model that the only scalar tetraquark systems to bind would be 
in the form of ,Of -molecules. This would predict just four states, the /o and a®, and not a 
complete nonet. The distiction between qq and qqqq states is most marked in the theorists' 
favourite world of large N c (where N c is the number of colours). Then the qq nonet in 
Fig. 4 becomes stable, while the tetraquark mesons merge into the continuum. Of course, 
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Fig. 4. The spectrum of scalar states given by a simple ideal nonet of qq and of qqqq mesons, where n = u,d, 
indicating which would be identified with the isoscalar <x and its isodoublet partner, the k, if the observed hadrons 
had these orthogonal compositions. 



we do not live in a world in which N c is large. Indeed, nowhere is this more apparent than 
in the scalar sector. Pelaez | 39 | 40 | jj as Q^^gd jjjg p- anc [ K as ^ j s varied and concluded 
that they are largely tetraquarks. With the addition of a glueball with mass co mputed on 
the lattice in a world without quarks to be in the region of 1.5-1.7 GeV |4JJ 42|43| we have 
our 19 "lightest" scalars. Of course, these states qq, "qqqq, (or even molecules) and glueball 
are not orthogonal to each other. They inevitably mix. So how can we try to distinguish 
which scalars have what composition? In particular, what is the composition of the lightest 
scalar H31 the <x? 



2. Correlation: photons probing hadron composition 

Two photon interactions can illuminate these issues. For orientation, let us presume we 
have extracted the two photon couplings of resonant states. Before looking at the enigmatic 
scalars, let us consider the two photon couplings of the tensor mesons, f% ai and fL that 
belong to an ideally mixed qq multiplet. Their radiative widths are proportional to the 
square of the average charge squared of their constituents, Fig. 5. The absolute scale of 
their couplings depends on dynamics: on how the qq pair form the hadron. In the non- 
relativistic limit, as with charmonia, this is simply related to the wavefunction at the origin. 
If we assume that these are equal for the tensor states, then we have the prediction that 

r(/2 -» 77) : T( fl2 -» yy) : F(/ 2 ' -» yy) = 25 : 9 : 2 . (1) 

Experiment Q i s j n reasonable agreement with this, given the uncertainties in extracting 
couplings from data (as we shall see). If we apply the same ideas to the pseudoscalars, 
n°, q, q' , though not ideally mixed, we would not be able to reproduce experiment^-^ 

F(7r° -> yy) : T(q -> yy) : T{rf -> yy) * 1 : 60 : 500 . (2) 



February 1,2008 8:4 WSPC/INSTRUCTION FILE 



review-mpla-rev 



Two photon coupling of the sigma 1 



Y 



2 



Y 




Fig. 5. Two photon decay rate of a quark model hadron is the modulus squared of the amplitude for yy to 
produce a qq pair and for these to bind by strong coupling dynamics. In present calculations these are treated with 
crude approximations. A genuine relativistic strong coupling approach is essential for the lightest pseudoscalars 
and scalars. 

While the non-relativistic quark model determines the intrinsic coupling to photons, it 
does not include the dynamics of binding the quark and antiquark into a meson, Fig. 5. 
A Lagrangian for the pseudoscalar-photon-photon interaction would introduce a factor of 
mass cubed for each meson. As noted by Hayne and Isgur 53 such factors are just what 
is needed to bring the quark model prediction into agreement with experiment of Eq. (2). 
However, what this teaches us is that we should not be using the non-relativistic quark 
model and correcting it by factors of hundreds for light quarks. Rather we need genuinely 
relativistic strong coupling calculations of such radiative widths. As we shall see, this need 
equally applies to the scalars. 

Now how do we extract such two photon couplings? Let us consider the exclusive pro- 
cess for which we have the most information, viz. yy — > tttt, with both charged and neutral 
pions, as displayed in Fig. 6. At low energy, the photon has long wavelength and interacts 
with the whole of the hadron. Consequently, it "sees" the charged pions, but not the neutral. 
7r°7r production is small. n + n~ is large (Fig. 6). How large is determined by the charge of 
the pion. However, as the energy increases (by just 1 GeV) the wavelength of the photon 
shortens sufficiently that it recognises that pions, whether charged or neutral, are made of 
the same charged constituents, namely quarks, and causes these to resonate. Since photons 
have spin-1, they readily produce the prominent tensor meson /2(1270), seen in Fig. 6. As 
the energy increases further, the photon probes not constituent quarks, but current quarks. 
Their interactions can be calculated perturbatively as shown by Brodsky and Lepage^! Al- 
lowing for a free normalisation, then these predictions for both the energy and c.m.angular 
distribution agree well with the latest data from Belle EH. What we are interested in here is 
the scalar signal underneath the large /2(1270) from nn threshold upwards, Fig. 6. 

To see how to extract this signal, let us focus on nn production close to threshold. 
Since pions are the Goldstone bosons of chiral symmetry breaking, their interactions at low 
energy can be computed perturbativelyE21 For yy — > n°7r°, the lowest order in chiral per- 
turbation theory (xPT) involves one loop gr aphs which predict a cross-section rising almost 
linearly with energy from threshold I 53 | 54 | as seen m pjg 7 This is in rather poor agree- 



ment with the only normalised data from Crystal Ball ^4 This was once a little worrying, 

feci 

since Maiani JJ described this process as providing a "gold-plated" prediction of^Pr, sug- 
gesting that the data might be wrong. However, one can calculate this cross-section almost 
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Fig. 6. Integrated c ross -section for yy — > nn as a function of cm. energy M(ror) from Mark II S3 Crystal 
Ball EZEH and CLEoEH. The n°n° results have been scaled to the same angular range as the charged data and by 
an isospin factor. Below are graphs describing the dominant dynamics in each kinematic region, as discussed in 
the text. 
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Fig. 7. Cross-section for yy — » n°n° integrated over | cos 8*\ < 0.8 as a function of the nn inva riant mass M(nn). 
The data are from Crystal BallEjSjThe line is the prediction of %PT at one loop (W)E2l The shaded band 
shows the dispersive prediction | 59 | 61 | 62 | described later — its width reflects the uncertainties in experimental 
knowledge of both Tin scattering and vector exchanges. 
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exactly at low energy, by noting that neutral pions can be produced by fi rst yy — » by 
the Born amplitude and then can interact to form the ttV final st ateES2EBE5]. That 
this cross-section is calculable from strong interaction dynamics makes it almost u nique 
among hadronic processes. This is because of Low's low energy theorem [ 60 | 59 | 6 1 1 6 2] ^his 
states that the amplitude for yn — > yn at threshold is exactly given by the one pion ex- 
change Born amplitude. Though this specifies the amplitude at just one kinematic point, 
the amplitude smoothly approaches this limit throughout the low energy region. This is be- 
cause of the closeness of the t and M-channel pion poles to the physical region for all three 
channels. This means that the Born amplitude modified by final state interactions controls 
the yy — > ttjt amplitude in the low energy region, where the pion pole is so much closer 
than any other exchange: the next nearest being the p and a>. Consequently, if we consider 
the analytic structure of the partial wave amplitudes, Fig. 8, they have a left hand cut gen- 
erated by crossed channel exchanges, in which essentially the part from s — to s =s -m 2 p 
is known from pion exchange, and a right hand cut controlled by direct channel dynamics. 
(As usual s is the square of the yy and nn cm. energy.) The discontinutity across this right 
hand cut is specified by unitarity: 

ImTiyy -> tttt; s) = ^ T(yy -> H n ; s)* T(H„ -> nn; s) , (3) 

n 

where T, T are partial wave amplitudes with definite isospin and the sum is over all al- 
lowed hadronic intermediate states H„, see the inset in Fig. 8. In the low energy region 

7171 
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Fig. 8. The complex i-plane structure of the yy — > Tin amplitudes, n labels the the left hand cut controlled 
by the pion exchange Born term, and p denotes where the vector exchanges contribute to the discontinuity. The 
discontinuity over the right hand cut is determined by the unitarity relation depicted. The nn intermediate state is 
effectively the only one possible from nn threhold up to the opening of the KK channel. As discussed in the text, 
this implies the yy — » nn amplitudes is "exactly" calculable in the region shown. 
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only nn is possible. Consequently, knowing the nn — > nn partial wave amplitudes from 
experiment, we can compute the yy — ♦ nn partial waves. 

To proceed with this calculation, let us first consider the S -wave yy — » nn amplitudes 
with isospin /, T'{s) — we will comment on higher partial waves later. Each of these 
amplitudes, with / = 0, 2, being complex has a phase (f>'(s) along the right hand cut, when 
s is above the two pion threshold, i.e. s > s t h = 4w^. In the elastic region, unitarity, 
Eq. (3) implies Watson's theorem, requiring the phase of each of these partial waves to be 
the same as the phase of the corresponding nn partial wave amplitude with the same spin 
and isospin (in principle, up to integer multiples of n. However, at nn threshold, the phases 
are equal to zero.) To implement this constraint we define the Omnes function, Q.'( s), which 
by construction has phase <p'(s): 

Q'(s) = exp 

Thus the yy — > nn S -wave amplitudes, T'is), can be written as P'(s) fl 7 (s), where P'(s) 
is a function which is real along the right hand cut with s > s t h- The phase, <p', is simply 
the phase-shift in the region of elastic unitarity, which applies up to KK threshold, since 
multi-pion channels are negligible below 1.2 GeV. Moreover, in the low energy region of 
interest where | s | ~ 0.25 GeV 2 , the unknown phase above 1 GeV affects the results little. 
This is checked by replacing the nn — > nn phase with that for nn — » KK. Such a change is 
equivalent to assuming the nn final state in the two photon process is only accessed through 
a KK intermediate state. Outside the narrow confines of the /o(980) region, this would be 
an extreme possibility. Nevertheless, the effect is small and included in the uncertainties we 
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Fig. 10. Cross-section for yy — > n + n integrated over |cos0*| < 0.6 as a function of the nn in variant mass 
M(nri). The data are from Mark IlEH. The shaded band shows the dispersive prediction | 59|6» 1 based on the 
Born term modified by calculable final state interactions in all partial waves. The width of the band reflects the 
uncertainties in experimental knowledge of both mz scattering and vector exchanges. 



quote. Representative input nn S -wave phases, <p'(s), for / = 0, 2 and the resulting Omnes 
functions are shown in Fig. 9. 

As already discussed the lightness of the pion means that Born amplitude not only 
equals the whole amplitude as s — > 0, and t, u — > m%, at the threshold for Compton scat- 
tering jti — > yn, but in the whole low energy region, as illustrated in Fig. 8. To implement 
this, let us begin by denoting the left hand cut contribution to the yy — > nn partial wave 
amplitudes generated by crossed channel exchanges collectively by £,'(s). As indicated in 
Fig. 8, pion exchange determines the discontinuity in the whole region > s > -Mh, 
beyond which other exchanges like p, a> start to contribute. While the Born term assumes 
pointlike couplings for the pion, any form-factor dependence only affects the left hand cut 
for s < -My, since it is vector masses, My, that set the scale for such charged radii. Con- 
sequently, the left hand cut from s — to s - -0.5 GeV 2 is precisely known and that is all 
we require to fix the amplitude in the region of the pole at s = sr = E 2 R in Fig. 2. 



To see how, let us construct the function G ! (s) = yF ! (s) - -C ! (s)j Q. l (sy l , which only 
has a right hand cut. Its discontinuity is £,'(s) sin0 / (s)/'|Q / (.s)|, which is accurately known 
at low energies. This information is embodied in a dispersion relation for the function G'(s) 
using a contour round its right hand cut and closed at infinity. While the behaviour of G'(s) 
means the integral at infinity converges with just one subtraction, it is more convenient for 
our purpose to ensure that the integrals are dominated by the known low energy regime of 
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< M 2 . This is achieved by making two subtractions: 

, , , s 2 , r°° , £V) sin^V) 

r'(s) = £ ! ( S ) + Curtis) + -q 7 (s) ds' y ; . (5) 



'Sift 

The constants c/ are specified by the QED low energy theorem and chiral dynamics. These 

two conditions apply to the amplitudes with the pions of specific charge (which are c ombi- 

B3 64 



nations of those with definite isospin). Low's theorem, as shown by Goldberger et au 
requires that the charged pion S -wave amplitude ( T c for yy — > n + 7T~), 



T c (s) -> B(s) + 0(s 2 ) as s -> 0, (6) 

where B(s) is the Born S -wave, while chiral dynamics requires that the neutral pion S -wave 
amplitude (T" for yy — > ttV) 

= at s = 0(m£) . (7) 

At one loop level in Chiral Perturbation Theory E31 T"(s) oc T{n + n~ — > 7r°7r°) and so 
places the Adler zero exactly at s — m 2 . at this order. However, its precise position hardly 
affects our results. An analogous procedure fixes t he higher waves too. Their threshold 
behaviour is enough to fix their subtraction constants^2|6JJ62(of these waves, the D-wave 
with helicity two is critical to the calculation of 

observablesfei£2l. 

In this way the yy — > nn cross-section in the low energy region is absolutely fixed. Pre- 
cision comes from the accurate determination of the nn S -wave amplitudes^|66| obtained 
by combining new results from decays like K e 4, J/tj/ — > (f>X and D — > nX (X = nn, KK) 
with the Roy equations, which fixes the pole at s — s R . This calculation reproduces the 
cross-section for the production of charged and neutral pions as measured by Mark II S3 
and Crystal Ball 321 respectively, in the low energy region with no free parameters, as dis- 
played in Figs. 7, 10. The bands shown delineate the uncertainties due to (i) different yy 
phases <p ! (s) above KK threshold and (ii) different positions of the Adler zero in Eq. (7). 
Notice that the cross-section is very nearly unique up to 450 MeV. That the yy — > nn 
amplitudes are reliably calculable at low energy from the Born amplitude modified by ex- 
perimentally determined final state interactions has two important applications. 



3. Radiation: the two photon coupling of the <x 

Of course, the / = nn phase and Omnes function, shown in Fig. 9, know about the cr-pole 
at s — sr deep in the complex plane close to both the right and left hand cuts of Figs. 2 
and 8. Not only can we determine the yy amplitudes T'{s) along the upper side of the right 
hand cut on the physical sheet where experiments are performed, but everywhere on this 
first sheet shown shaded in Fig. 2. In particular, we can determine the 1 = amplitude at 
s = sr, marked by the dot and cross in Fig. 2. 

The right hand cut structure, Fig. 8, of the yy — > nn amplitude mirrors that of the 
corresponding hadronic amplitude, T 1 , for nn — > nn in the region of elastic unitarity, so 
that: 

T'is) = a '(s)T'(s) . (8) 
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This satisfies Eq. (3) if a'(s) is a real function for s > 0. The function, a'(s), represents the 
intrinsic coupling of yy — > nn, while the nn — > nn amplitude, T , describes the final state 
interactions, which colour and shape the electromagnetic process 21 The function a'(s) is 
determined by crossed-channel dynamics — the N of N/D for the cogniscienti. 

At s = sr on the first sheet, the amplitude T I=0 (s) = i/2p(s), since the S -matrix 
element vanishes at this point. p(s) is, as usual, the phase-space factor p(s) = y/l - s,i, /s. 
The dispersion relation on the first sheet then determines the coupling function ch(sr), which 
not having a right hand cut, has the same value on the second sheet of Fig. 2. 

Let us introduce subscripts to label the sheets / and //, while the superscripts continue 
to denote isospin. In the neighbourhood of the pole on the second sheet, the yy — > nn 
S -wave amplitude is given by 

2 

Tf t {s) - , while Tfjis) - . (9) 

SR - S Sr - S 

Even if factorised residues are not strictly appropriate for such a very short-lived state, 
Eqs. (10,11) below provide a physically meaningful and unambiguously defined EH yy 
width. a (sr) determines the ratio of g 7 lg„ for the isoscalar resonance. Now the hadronic 
amplitude on sheet I is related to that on sheet II by 

' ' +2ip , (10) 



T n (s) Ti(s) 
so that 



2 (s- S r)F»(s) 2 
?y = llm — ~ ■ (11) 



Ms 



(7» - i/2p) 



The 1 = amplitude T° encodes the key nn final state interactions a nd so a utomatically 
incorporates the effect of the cr. Indeed, the representation, cited above ^ for 7"~° , has a 
pole at E = 441 — i 227 MeV, within the error ellipse found by Caprini et al. 021. Inputting 
this hadronic amplitude on sheet I into the present dispersive calculation then determines 
the residue of the cr-pole in yy — > nn, and in turn specifies EES its radiative width to be 

a 2 \p(s R ) gl\ 

Y(o-^yy) = y - = (4.1 + 0.3) keV , (12) 

as found in Ref. 69. 

The low energy theorem and knowledge of nn scattering has allowed us to determine 
the low energy yy — > nn partial waves that accurately sum to describe both the charged and 
neutral data of Figs. 6, 7 and 10. In these data the broad scalar is far from apparent. The 
near threshold peak in the charged cross-section is due to one pion exchange Born term, 
modified by final state interactions. The presence of this Born term, means there are large 
7 = and 1 = 2 components in the yy amplitude. These essentially add in the charged 
pion channel and subtract in the neutral one. It is this that makes the "cr" component far 
from obvious in the data. Nevertheless, one naturally asks: how does the large radiative 
width for the cr of 4 keV square with the size of the yy — > nn cross-sections shown in 
Figs. 7, 10. Simplistic reasoning would totally ignore the key requirement that final state 



I 



I 
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I i i r 




C 



b i i i d 

Fig. 11. yy — » nn S -wave amplitudes, F, at 400 MeV with definite isospin and with definite charges as indicated 
by the superscripts. B is the Born amplitude S-wave for comparison. OC and ON define the directions of the 
charged and neutral pion amplitudes as given by the vector sums described in the text. 



interactions shape the nn distribution in a well-defined way. Then one would say that, with 
no Born contribution in the neutral pion channel, the cross-section should directly reflect 
the appearance of resonant structures. If this is the cr, then one can read off, from the 
observed cross-section in Fig. 7 of 10 to 12 nanobarns, a yy width an order of magnitude 
smaller than we have deduced. However, this is too naive. 

The reason is that final-state interactions are crucial. These apply quite differently to 
the 7 = and 2 yy amplitudes. In hadronic reactions 7 = 2 amplitudes, being exotic in the 
quark model, are much smaller than those with / = 0. In contrast in the two photon process 
both 1-0,2 are equally important. The <x appears in the 7 = amplitude, and this can 
only be separated from data by analysing yy — > n + n~ and 7r°7r° tog 

As we have seen what is really happening is that the Born amplitude is modified by 
final state interactions to ensure Watson's theorem is satisfied. As a result the 7 = and 2 
S -wave amplitudes are no longer real and exactly cancelling in the neutral channel. A 
vectorial representation of this is shown at 400 MeV in Fig. 1 1 . The 7 = component 
has the phase of 7 = S -wave nn scattering, while that with 7 = 2 has the phase of the 
corresponding isotensor S-wave. In Fig. 1 1 the vector OC (T +T 2 1 V2) is V3/2 times the 
charged channel S -wave, F°, while the vector ON (T - yflT 2 ) is - V3 times the neutral 
one, F". One sees that the square of the neutral channel S -wave (which dominates its cross- 
section) is a factor of 12 smaller than the modulus squared of the pure 7 = S -wave. It 
is in this (bigger) amplitude that the cr is to be found. This delivers an 7 = cross-section 
averaged across the cr consistent with a 4 keV width determined from the pole residue. 

There are, of course, other calculations of the two photon width of the cr in the litera- 
ture LZQ] a particularly recent one is that of Fil'kov and Kashevarovl^l] These authors fit 
the data in Fig. 6 with amplitudes that satisfy dispersion relations, with a resonant shape 
for the imaginary part of the cr-component. In Ref. 72 they have three fits, which deter- 
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T(yy) keV 


composition 


predictions 


author(s) 


(uu + dd)/^/2 


4.0 


Babcock & Rosner^S 


Is 


0.2 


Barnes ^31 


[ns][ns], n = (u,d) 


0.27 


l7Sl 

Achasov et al. 11 ^ 


KK 


0.6 
0.22 


Barnes 
Hanhart et al. 



Table 1 : Radiative widths of scalars in different modellings of their composition. 



mine Breit-Wigner parameters for the <x spanning a whole range of values from 449 to 
616 MeV for the mass and 542 to 1655 MeV for the total width. These give corresponding 
two photon widths from 0.26 to 0.99 keV, with 0.62 keV favoured in their later Ref. 71 — 
much smaller than our result of Eq. (12). However, the treatment of Fil'kov and Kashe- 
varov does not respect the key unitarity relation of Eq. (3). As just discussed in connection 
with Figs. 7, 10, 11, ignoring this leads to erroneous conclusions. The way the phase of the 
7 = 7 = component changes in the complex energy plane, Fig. 3, is crucial in the linking 
the pole parameters to the behaviour of the relevant amplitudes on the real energy axis from 
E = to E = 0.6 GeV, in particular. 

What does our radiati ve width of (4.1 + 0.3) keV mean for the composition of the cr? 
Simple prediction iP'^^'^^"^ are shown in Table 1. We see that the result in Eq. (12) 
seems to affirm a simple uu , dd composition. However, the calculation of the qq composi- 
tion is not really robust. It can be checked by a relation provided by the naive quark model 
for states in the same L-band. Here, tensors and scalars of the same quark constitution, with 
the same spin (5 = 1) and the same orbital angular momentum (L = 1), satisfy GB|79|. 

F(0 ++ -> yy) = 15 lmo\" 
T(2++ yy) 4 \m 2 j ' 

where the exponent n in the mass term depends on the relevant binding shape of the in- 
terquark potential: n = 3 for the short-distance Coulombic component, and n — > for the 
linearly confining part. With the <r and fx masses differing by more than a factor of 2, the 

or) 

result is very senstive to this exponent n. Moreover, Li et al. I2iil have estimated that this 
relation can be changed by up to a factor of 2 from relativistic corrections. 

IO j I 

Alternatively, Narison l2JJ has proposed that the cr is largely a glueball. Perturbative 
arguments would suggest a small two photon coupling for such a state. This is certainly 
difficult to reconcile with the calculated rate of 4 keV. Once again we send out a plea for 
a genuine strong coupling calculation of the scalar radiative widths with different com- 



I 



I 
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positions qq, qqqq or glueball to compare with the result of Eq. (12) and the coupling in 
Eq. (11), from which it is derived. 

An important check of the calculation presented here and displayed in Figs. 7,10 would 
be a more precise measurement of the charged and neutral pion differential cross-sections 
at low energy. This is the one part of the proposed DAcDNE physics programme that has 
not yet been implemented. While the luminosity is potentially higher than either PETRA or 
PEP, where the Mark II and Crystal Ball data were tak en, the c onfiguration of the DAONE 
machine requires the on-going electrons to be tagged 1 1 1 I This is a feature that should 
be taken up in any upgraded DAONE-II. At higher energies, i.e above 800 MeV, data are 
now becoming available from Belle, with the charged cross-section already published^] 
and neutral data in the pipeline. These will provide results, both integrated and differential, 
with high statistical precision that will impact on the first accurate determination of the yy 
coupling of the higher mass scalars, the /o(980), /o(1370) and /o( 15 10). Detailed amplitude 
analyses will be required to separate out the / = J — component of these cross-sections 
cleanly, as discussed recently in Ref. 86. In addition Belle will have data on the yy — > jfiq 
channel, which will be key to unravelling the composition of the ao(980) and ao(1430) and 
so learning about which of the isoscalars is their companion. Which of these, beyond the 
<x, make up the Higgs sector of QCD is not yet clear. 

Modelled on the way the uu , dd condensates break chiral symmetry, heavy flavour con- 
densates may similarly break electro-weak symmetry — see for instances Refs. 87-89. The 
first priority is to determine the location of the corresponding scalar(s). After that the cor- 
relation of couplings to the many decay channels must be found, Two photon couplings so 
essential for the study of the Higgs sector of chromodynamics, are likely to prove just as 
important for illuminating the Higgs of flavourdynamics at future colliders. The radiation 
of photons is the key to solving the riddle of the scalars. 
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